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ABSTRACT 

We report the discovery of three binary millisecond pulsars during the Parkes Multibeam Pulsar 
Survey of the Galactic Plane. The objects are highly recycled and are in orbits of many tens of 
days about low-mass white-dwarf companions. The eccentricity of one object, PSR J1853+1303, is 
more than an order of magnitude lower than predicted by the theory of convective fluctuations during 
tidal circularization. We demonstrate that, under the assumption that the systems are randomly 
oriented, current theoretical models of the core-mass-orbital-period relation for the progenitors of these 
systems likely overestimate the white-dwarf masses, strengthening previous concerns about the match 
of these models to the data. The new objects allow us to update the limits on violation of relativistic 
equivalence principles to 95% confidence upper limits of 5.6 x 10 -3 for the Strong Equivalence Principle 
parameter |A| and 4.0 x 10 -20 for the Lorentz-invariance/momentum-conservation parameter |qt3 | . 

Subject headings: pulsars: individual: PSR J1751-2857, PSR J1853+1303, PSR J1910+1256 - 
stars:binaries — relativity 



1. INTRODUCTION 

Millisecond radio pulsars are the product of an 
extended period of mass and angular momentum 
transfer to a neutron star (NS) from an evolv- 
ing companion star. This "recycling" s c enario 

was proposed iBisnovatvi-Kogan fe Komberd 119741 
ISmarr fc Blandfordl 1197ft lAlmr et a,U 11 9821) shortly 
after the discovery of t h e firs t binary and r ecycle d 
pulsars tHu lse fc Tavlorl 119751: IB acker et alJ 11982ft . 
Overviews of the mass transfer process are given in sev- 
eral places (e.g., iBhattacharva fc van den HeuvellT991l: 



Phinnev fc Kulkarni 1994; ITauris fc van den Heuvell 
It has become clear that there are in fact 
several sub-classes of recycled pulsars, the most obvious 
distinction being between those that have NS versus 
white-dwarf (WD) companions. Even within the latter 
category and focusing on Galactic field binaries only, 
a wide range of companion masses and evolution- 
ary histories are represented. One sub -group is the 
"inte r mediate-mass binary p ulsars" (e.g., iCamilo et alJ 
l200lt lEdwards fc Bailed l200l|) which have spin periods 
of tens of milliseconds and/or companions that are likely 
massive (CO or ONeMg) WDs. Many of these systems 
may have experienced an ultra-high mass transfer rate, 
or else undergone a period of common-envelope evolution 
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llvan den Heuvell 119941: iTaam et alJ 120001 ITauris et all 
|2000(k Among pulsars with lower-mass Helium WD 
companions, there is another split. Systems whose 
initial orbital periods were le ss than the "bifurcation" 
period of about 1 or 2 days ijPvlvser fc Savoniie] 119881: 
lErgma et alJ 11998^) will see their orbits shrink through 
magnetic braking and gravitational radiation, ultimately 
forming a low- mass binary pulsar in a tight orbit. 
Systems with longer initial orbital periods undergo 
stable, long-lived mass transfer via an accretion disk, 
once the companion star has evolved onto the giant 
branch. The resulting systems have long orbital periods 
of several days or more. The prototype of this "wide- 
orbit binary millisecond pulsar" (WBMSP) group, 
PSR B1953 +29, was one of the first recycled pulsars 
discovered (jBoriakoff et alJ I1983J) : over the years the 
number of these pulsars with orbital periods greater 
than 4 days has grown to 18. 

The WBMSPs are the best-understood class of pulsar- 
WD binaries. For instance, the companion mass and 
orbital separation are thought to follow the "core-mass- 
orbital-period" (Pb—mi) relation, in which the mass of 
the core which eventually forms the white dwarf is di- 
rectly related to the size of the envelope of the Roche- 
Lobe-filling giant star and hence the orbital radi us (e.g., 
Rapp anort et all 11995: Ta uris fc Savoniie] 119991) . Mea- 
surements of companion masses to date (via pulsar tim- 
ing or optical spect roscopy) indicate that this relation is 
fairly well satisfied Kj^i j^jdjjl^^ aTl 



I200"ll ISr^Ter et aljISoOft Ivan Kerkwiik et alJEOOflT aT 

though we shall discuss this issue further below. Fur- 
thermore, although the orbit becomes tidally circularized 
during the companion's giant phase, fluctuations of con- 
vective cells in the gian t's envelope for ce the eccentric- 
ity to a non-zero value l)Phinne vl ll992f) : thus, to within 
an order of magnitude or so, the eccentricities of these 
WBMSP systems can be predicted from the orbital pe- 
riods. 

More examples of WBMSP pulsars are needed 
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to further test these predictions, and also to pro- 
vide additional constraints for population synthe- 
sis efforts (e.g . . iPqrteeies Zwart fc Yungelsonl Il998t 
IWillems fc Kolhl 120021 12003 IPfahl et all 12003ft . These 
objects are also valuable for constraining departures from 
general relativity (G R) in the form of equival ence prin- 
ciple violations (e.g.. iDamour fc Schaferlll99l . Finally, 
some pulsar- WD binaries permit measurement of NS and 
WD masses through relativistic or geometric timing ef- 
fects; thus new systems potentially add to the pool of 
objects that can be used to constrain theories of the NS 
interior and/or the amount of mat ter transferr ed in dif- 
ferent evolutionary processes (e.g.. IStairsll2004ft . In this 
paper we report the discovery of three more WBMSPs 
during the Parkes Multibeam Pulsar Survey. In § [21 wc 
describe the search and follow-up timing observations. 
In § 13.11 wc describe the characteristics of the new pul- 
sars and relate them to evolutionary theory including a 
comparison with the predictions of the P^-m^ relation. 
In S 13. 21 we use the ensemble of pulsars with white-dwarf 
companions to set new, stringent limits on equivalence 
principle violations. Finally, in § 0] we summarize our 
results and look to the future. 

2. OBSERVATIONS AND DATA ANALYSIS 

The Parkes Multi beam Pulsar Survey (e.g., 
iManchester et all 120011) used a 13-beam receiver on 
the 64-m Parkes telescope to search the Galactic Plane 
(|6| < 5°, 260° < I < 50°) for young and recycled 
pulsars. Observations were carried out at 1374 MHz 
with 96 channels across a 288 MHz bandpass, using 
35-minute integrations and 0.25 ms samp ling. More 
than 700 pulsars have been discovered (e.g. JHobbs et al.1 
12004ft , nearly doubling the previously known population. 
The survey processing includes dedispersion of the data 
at numerous trial dispersion measures, followed by a 
periodicity search u sing Fast Fourier Tran sforms and 
harmonic summing (Ma nchester et al.ll2001ft . Recently, 
the entire data volume has been reprocessed using "ac- 
celeration" searches for pulsars in fast binary orbits, as 
well as improved interference excision techniques. This 
has resulted in a number of new binary and millisecond 
pulsars (Faulkner et al. 2004, 2005), including two of 
the objects described in this paper. 

PSR J1751-2857 has been observed since MJD 51972 
with Parkes at 1390 MHz using a 512-channel filterbank 
across 256 MHz bandwidth and 0.25 ms sampling and 
with the 76 m Lovell Telescope at Jodrell Bank Obser- 
vatory at 1396 MHz using a 64-channel filterbank across 
64 MHz with 0.13 ms sampling. 

PSRs J1853+1303 and J1910+1256 were discovered in 
observations taken on MJDs 52321 and 52322, respec- 
tively, and identified as candidat es using the REAP ER pro- 
gram selection procedure (iFaulkner et aljr2004ft . Both 
were confirmed as pulsars on MJD 52601. Subsequently, 
they have been timed at Parkes and Jodrell Bank us- 
ing the systems described above. Both have also been 
regularly observed using the 305 m Arecibo telescope 
in Puerto Rico. These observations initially used one 
100 MHz Wideband Arecibo Pulsar Processor (WAPP) 
in a fast-sampled "search" mode, and, since 2004 Feb., 3 
WAPPs in an online folding mode, in all cases with 256 
lags and 64 /isec sampling. The single- WAPP observa- 
tions were centered at 1400 MHz and when 3 WAPPs 
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Fig. 1. — Timing residuals and full-period profiles at 1400 MHz 
for each of the three new pulsars. Top: PSR J1751-2857. Middle: 
PSR J1853+1303. Bottom: PSR J1910+1256. The profile for 
PSR J1751— 2857 was obtained with Parkes; the other two with 
Arecibo. For each pulsar, the TOAs with large scatter are from 
Jodrell Bank for all pulsars and from Arecibo for PSR J1853+1303, 
and represent eras when the timing solutions were poorly known. 



were used they were centered at frequencies of 1170, 
1370 and 1470 MHz. Some of the single- WAPP obser- 
vations display timing systematics similar in both pul- 
sars relative to the contemporaneous Parkes data; as 
this likely reflects instrumental errors in the fairly new 
WAPP, these data points have been left out of the tim- 
ing analysis. Some of the Arecibo observations were flux 
calibrated using a pulsed noise diode of known strength; 
the resulting calibrated profiles are used to determine 
the flux densities for PSRs J1853+1303 and J1910+1256. 
The flux density for PSR J175 1-2857 was determ ined 
using the procedure outlined in lHobbs et alJ (|2004ft . 

The data from each telescope were dedispersed and 
folded modulo the predicted topocentric pulse period; 
this was accomplished off-line for the Parkes data and 
Arecibo search-mode data and on-line for the Jodrell 
Bank and Arecibo folding-mode data. A Time-of- Arrival 
(TOA) was determined for each observation by cross- 
correlation w ith a high signal-to-noise standard template 
l)TavlorJll992ft . The timestamp for each observation was 
based on the observatory time standard, corrected by 
GPS to Universal Coordinated Time (UTC). The JPL 
DE200 ephemeris (|Standisblll990() was used for barycen- 
tric corrections. The timing solutions were found using 
the standard pulsar timing program tempo 8 , with uncer- 
tainties containing a small telescope-dependent amount 
added in quadrature and scaled to ensure xt — 1- 
The resulting pulsar parameters are shown in Table ^ 
while the timing residuals and standard pulse profiles 
are shown in Figure ^ We note that the correctly folded 
Arecibo data have rms residuals on the order of 1 /isec 
per WAPP for PSR J1910+1256 and 1-2 /isec per WAPP 
for PSR J1853+1303 for roughly 30-minute integrations. 

3. DISCUSSION 
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TABLE 1 

Parameters for the new pulsars 



Parameter J1751-2857 

Right ascension, a (J2000) 17 h 51 m 32! 6965(2) 

Declination, 5 (J2000) -28° 57' 46'.'50(3) 

Pulse period, P (ms) 3.9148731963690(6) 

Period derivative, P (ss^ 1 ) 1.126(4) X 10~ 20 

Epoch (MJD) 52560.0 

Dispersion measure (pccm" 3 ) 42.808(20) 

Orbital period, P h (d) 110.7464576(10) 

Projected semi-major axis, x (lt-s) 32.528221(9) 

Eccentricity, e 0.0001283(5) 

Longitude of periastron, ui (deg) 45.52(19) a 

Epoch of periastron, To (MJD) 52491. 58(6) a 

Data span (MJD) 51808-53312 

Number of TOAs 168 

Weighted RMS timing residual (/is) 28.5 

Flux Density at 1400 MHz, S1400 (mjy) 0.06(2) 

Derived Parameters 

Galactic longitude, I (deg) 0.65 

Galactic latitude, b (deg) —1.12 

Distance (Cordes fc Lazio 20021 (kpc) 1.4 

Mass function, f\ (M Q ) 0.003013034(2) 

Minimum companion mass, m,2 (Mg) b 0.19 

Surface magnetic field, B = 3.2 X 10 19 \J PP (G) 2.1 X 10 8 

Characteristic age, t c = P/2P (Gyr) 5.5 



J1853+1303 



J1910+1256 



18 h 53 m 57^31827(8) 
13° 03' 44('0884(17) 
4.0917973806819(14) 
8.85(10) xlO- 21 
52972.0 
30.5702(12) 
115.6537868(4) 
40.7695200(10) 
0.00002369(9) 
346.63(9) a 
52890.25(3) a 
52606-53337 
140 
2.9 
0.4(2) 



19 h 10 m 09! 70041 (6) 
12° 56' 25^5276(6) 
4.9835839397055(12) 
9.77(7) xl0~ 21 
52970.0 
38.0650(7) 
58.46674201(9) 
21.1291045(6) 
0.00023022(6) 
106.001(ll) a 
52968.4474(18) a 
52602-53337 
183 
1.8 
0.5(1) 



44.87 
5.37 
1.6 

0.0054396358(4) 
0.24 
1.9 x 10 s 
7.3 



46.56 
1.80 
1.9 

0.0029628402(2) 
0.22 
2.2 X 10 8 
8.1 



Note. — Figures in parentheses are uncertainties in the last digits quoted, which are twice the formal errors reported 
by TEMPO after scaling the TOA uncertainties to obtain a reduced-^ 2 of about 1.0. 

a The parameters lu and To are highly covariant. Observers should use for J1751— 2857: lu = 45.523832°, To = 



are 

52491.578696239; for J1853+1303: 
= 52968.447431428. 
b Assuming a pulsar mass of 1.35 Mq. 



346.630447°, T = 52890.248760182; for J1910+1256: 



106.001079°, T 




10 100 

P, (days) 



Fig. 2. — Eccentricity e vs. orbital period Pb for the pulsar- WD 
systems thought to be described by stable mass transfer and the 
Pb— m2 relation. The three new pulsars are labeled and indicated 
by triangles. The s olid and dashed curves illustrate the eccentricity 
ranges predicted by Phinncv 1 1992) as a function of orbital period; 
95% of pulsars should fall within this range and 90% of the observed 
systems do. The dotted line has P 2 oc e, indicating the figure-of- 
merit for tests of the Strong Equivalence Principle (SEP). 



3.1. Pulsar Parameters and Evolution 



Each of the new pulsars is in an orbit of several tens 
of days with a companion of minimum mass of 0.2 to 
0.25 Mq. The eccentricities of PSRs J 175 1- 2857 and 
J1910 +1256 agree well with the predictions of IPhinnevI 
(1992); however, that of PSR J1853+1303 is lower than 
the predictions by more than an order of magnitude 
(Fig- 13 ■ As it has a low-mass companion and its spin 
period indicates that it is highly recycled, there is little 
other reason to believe that its evolution proceeded in 
any unusual fashion, so its low eccentricity may simply 
reflect the natural scatter in the population. 

There are now 21 pulsars with low-mass WD compan- 
ions whose orbital characteristics should be determined 
by the P^-m.2 relation (Tabled- To date, however, there 
are only a handful of observational tests of this relation: 
three Shapiro-delay timing measurements and one opti- 
cal WD spectr um for systems with o r bital periods greater 
than 2 days ijSplaver et al.l 120051: Ivan Kerkwiik et all 
12005(1 . In order to judge the agreement of the whole 
population with the theory, therefore, we need to use 
statistical argu ments based on the observed ma ss func- 
tions. We follow rThorsett fc ChakrabartyJ l(1999[) in con- 
sidering only those pulsars with orbital periods greater 
than about 4 days, excluding those in the 2-4 day range 
as being too close to the limits of applicability of the 
relation. With this many pulsars, it becomes possi- 
ble to examine whether different subgroups are equally 
well described by the relation. To do this, w e adopt 
the approach of lThorsett fc Chakrabartvl l)1999D (Fig. 4 
and related discussion) by assuming a range of pulsar 
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TABLE 2 

The 21 pulsars with low-mass WD companions and orbital periods of more than 4 DAYS 



T— 

Pulsar 


— — — - — 

P (ms) 


P b (days) 


e 


u>{ ) 


T") TV IT T~>A 

FM HA 


Tit, f T"\ _ _ 

FM Dec 




h 


yj , c 

K.eierCIlCCS 






(mas/yr) 


(mas/yr) 


( M e) 




J0407+1607 


25.702 


669.0704 


0.0009368(6) 


192.74(2) 






0.002893 


1 


J0437-4715 


5.757 


5.7410 


0.0000191686(5) 


1.20(5) 


121.438(6) 


-71.438(7) 


0.001243 


2,3 


T1 045— 4509 


7.474 


4.0835 


00001 97(1 3") 


243(4) 




6(1) 


001 765 


4 5 


71 455 — 3330 


7.987 


76.1746 


0001 697(3) 


223 8(1 ) 


5(6) 


24(12) 


006272 


6,5 


J1640+2224 


3.163 


175.4606 


0.000797262(14) 


50.7308(10) 


1.66(12) 


— 11.3(2) 


0.005907 


7.8 


J1643— 1224 


4.622 


147.0174 


0.0005058(1) 


321.81(1) 


3(1) 


-8(5) 


0.000783 


6.5 


T1 709+231 3 


4.631 


22.7119 


00001 87(2) 


24.3(6) 


3 2(7) 


9 7(9) 


007438 


7.9 


J1713+0747 


4.570 


67.8255 


0.0000749406(13) 


176.1915(10) 


4.917(4) 


-3.933(10) 


0.007896 


10,11 


J1732— 5049 


5.313 


5.2630 


0.0000098(20) 


287(12) 






0.002449 


12 


J1751-2857 


3.915 


110.7465 


0.0001283(5) 


45.52(19) 


_ 


_ 


0.003013 


This work 


J1804-2717 


9.343 


11.1287 


0.000035(3) 


160(4) 






0.003347 


13 


J1853+1303 


4.092 


115.6538 


0.00002369(9) 


346.63(8) 






0.005440 


This work 


B1855+09 


5.362 


12.3272 


0.00002170(3) 


276.39(4) 


-2.899(13) 


-5.45(2) 


0.005557 


14,15 


J1910+1256 


4.984 


58.4667 


0.00023022(6) 


106.001(11) 






0.002963 


This work 


J1918-0642 


7.646 


10.9132 


0.000022(4) 


234(11) 






0.005249 


12 


B1953+29 


6.133 


117.3491 


0.0003303(1) 


29.55(2) 


-1.0(3) 


-3.7(3) 


0.002417 


16,17 


J2016+1948 


64.940 


635.039 


0.00128(16) 


90(5) 






0.009112 


18 


J2019+2425 


3.935 


76.5116 


0.00011109(4) 


159.03(2) 


-9.41(12) 


-20.60(15) 


0.010687 


19,20 


J2033+1734 


5.949 


56.3078 


0.00012876(6) 


78.23(3) 


-5.94(17) 


-11.0(3) 


0.002776 


21,15 


J2129-5721 


3.726 


6.6255 


0.0000068(22) 


178(12) 


7(2) 


-4(3) 


0.001049 


13,5 


J2229+2643 


2.978 


93.0159 


0.0002556(2) 


14.42(5) 


1(4) 


-17(4) 


0.000839 


22,17 



References. — 1. ILorimer et"a!1 120051 ) 2. I.Iohnston et alJ ll993l). 3 Ivan Straten et alJ 120011) . 4. IBailes et ail 119941) 
5. IToscano et aU 11993) ■ 6. ILorimer et alj 119951). 7 . Foster et all 119951). 8. I Lohmer eTaLI 120051). 9. ILewandowski et alJ 
120041). 10. Foster et"al1 119931). 11 lijplaver etaT] 120051). 12. [Edwards fc Bailesi 120011). 13. ILorimer et alj 119961). 14 
ISegelstein et alJ 11986))~15. ISplaverl 120041) 16. IBoriakoff et alJ 119831). 17. IWolszczan et alj 120001) . 18. INavarro erakl 
12003) . 19. INice et alj 119§3T) . 20. INice et alJ 120011) . 21. IRav et alJ JT996I) . 22. ICamilo et alJ 119961) 




Fig. 3. — Cumulative probability distributions for the mea- 
sured mass functions p(fi < /) (solid lines) and the median 
predicted values of cos i ( dashed lines) for the 21 binary sys- 
tems listed in Table after IThorsett fc Chakrabartvl 119991) . Fig. 
4. In each panel, the thinnest lines incorporate all systems, the 
medium-weight lines only those systems with J\> < 50 days, and 
the heaviest lines only those systems with > 50 days. Panel 
(a) ass umes m2 is drawn unif ormly from the range predicted 
by the Tauris & Savoniic (1999) Pb— rri2 relation and that mi is 
drawn fr om a Gaussian distributi on cent ered on 1.35MQwith width 
0.04M^ fThorsett &; Chakrabarty|19991) . Panel (b) draws mi from 
the same ra nge, but assumes mi i s drawn from the range pre- 
dicted by the RaDpaport ct al. ( 1995) relation. Panel (c) as (a), but 
drawing m\ from a Gaussian distribution centered on 1.75M0with 
wid th 0.04Mp,. Pane l (d) a s (b) , but with m-2 limits given by 
the ITauris & Savo niic ( 1999) fits to the [Rappa port et a U 119951) 
results. The straight lines indicate the cumulative probability for 
a uniform distribution. 



masses (e.g., 1.35±0.04M Q , which was a good match 
to the set of pulsar masses m easured at the time of 
IThorsett fc Chakrabartvl <ll999|)'l and the (uniform) P h - 
m ? relation specified by ei ther ITauris fc Savoniid l|1999fl 
or iRappaport et"aT] l)1995j) . For the first test, for each 
pulsar we assume a uniform distribution in cos i and sim- 
ulate a large number of systems, finding the probability 
p(fi < /) that the simulated mass function / is above 
the observed value f\. For the second test, we assume the 
observed /i and find the median predicted value of cosi, 
again simulating a large number of systems. The cu- 
mulative probability distributions for both cosi (dashed 
lines) and p(fi < f) (solid lines) are displayed in Fig. 
for all 21 pulsars, and for those with orbital periods less 
than (8 systems) and greater than (13 systems) 50 days. 
For comparison, the straight lines indicate the cumula- 
tive probability for a unifo rm distribution. 

The m,2 limits from the ITauris fc Savoniiel l)1999fl re- 
lation (panel (a) of Fig. |2J are given by the Pop. I 
and Pop. II fits to their simulation results; they find 
a spread in Pb of a factor of about 1.4 aroun d their 
median value for a given m 2 . Ra ppaport et al.1 l)1995fi 
find a spread of about 2.4 in Pb for any given TO2, but 
consider this to cover roughly the full range of possi- 
ble values; we ass ume a range of \/2A ^ 1.6 will be 
comparable to the ITauris fc Savoniiel l)1999) ranges and 
show the corresponding results in panel (b) of Fig. 
ITauris fc Savoniiel (119991) a lso provide their own fits to 
the lRappaport et al.1 1)1995) simulations, and we evaluate 
these fits in panel (d) o f Fig. |31 It is important to note 
that these lTauris fc jj ji avonije| <|199ffi fits do not cover the 
full spread of the Rappapo rt et alJ l)1995t) orbital peri- 
ods, favoring the lower periods at any given mi- Thus it 
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is perhaps not surprising that panel (d) indicates higher 
mass estimates in g eneral than panel (b). 

We find that the iTauris fc Savoniid {1999ft Pb~m 2 re- 
lation is incompatible at the 99.5% level (according to 
a KS-test) with a uniform distribution of cos i if the 
pulsar masses are drawn from a Gaussian distribution 
centered on 1.35M Q with width 0.04M Q . Better agree- 
ment with uniformity in cosi (at the 50% level) can 
be reached if the pulsar masses are very large on av- 
erage (e.g., 1.75±O.O4M ; panel (c) of Fig. EJ, a sit- 
uation not supported by observational evidence. The 
iRannanort et alJ l|1995ft relation appears to be in slightly 
better agreement with uniformity in cosi, though it is 
clear from Fig. [3] that this occurs because of a tendency 
to underestimate the companion masses for short-period 
systems and overestima te those for long-period systems . 
We note that, although iThorsett fc Chakrabartvl (|1999f> 
favor using the iRappaport fc Jossl l)1997fl version of the 
Pb-m 2 relation for m 2 < 0.25 M Q , this predicts ex- 
tremely low masses (~ 0.10 Mq) for the shortest- Pb 
systems, which is in conflict with the observed com- 
panion mass of 0.236 ± 0.017 M s for PSR J0437-4715 
l)van Straten et all 120011) . Using this revised relation 
would lower the estimates of the companion masses in the 
short-Pb systems even further; this in combination with 
the higher estimates for the long-Pb systems appears to 
have been responsible for t he ove rall good agreement that 
IThorsett fc Chakrabartvl l)1999l) found for the combined 
relation with NS masses of 1.35±0.04M Q and a uniform 
distribution in cosi. Thus, assuming that the cosines of 
the system inclination angles are in fact uniformly dis- 
tributed and that most NS masses are near 1.35 Mq, it 
appears that the existing forms of the Pb—rri2 relation 
tend to overestimate companion masses for long-period 
systems, while providing conflicting results for the short- 
period systems. 

The tendency to underesti mate masses in the short - 
Pb case was in fact noted bv IRappaport et all {l995), 
although they included in their analysis systems now 
considered to be "intermediate-mass" binaries (such as 
PSR J2145— 0750) havi ng different evolutionary histo- 
ries. ITauris! <jl996) and ITauris fc Savoniid Jl999) com- 
ment on the poor match of the (then 5) known WBM- 
SPs to the higher theoretical predictions of m-i . With the 
larger number of systems now k nown, the conclusion of a 
poor match seems inescapable. Rappapo rt et alJ (1995) 
note that while the relationship between core mass and 
luminosity for red giants is well understood, the relation 
between mass or luminosity an d radius is looser (see also 
IThorsett fc Chakrabartvll999r i. with uncertainties in the 
companion's initial chemical composition and the convec- 
tive mixing-length parameter; this may explain our re- 
sults. It appears more theoretical work will be required 
to derive models that better match the data. 

The only long-orbital-period system wit h a timing test 
of the P h -m 2 relation is PSR J1713+0747 ijSplaver et all 
2005), and the meas ured companion mass is in fact 
slightly lower than the ITauris fc Savoniiel 1)199 9) predic- 
tion. Mass measurements or constraints in more sys- 
tems, by timing or by optical spectroscopy of the WD 
companions, will be needed to confirm or refute our 
present conclusions. The new systems will likely lend 
themselves to observations of geometrical effects such as 
the change in apparent semi-major axis due to motion of 



the pulsar and/or the Earth l|Kopeikinlll995l Il996ft . and 
J1910+1256 in particular has sufficient timing precision 
that it may be possible to measure Shapiro delay in this 
system. 

3.2. Equivalence Principle Violations 

The WBMSPs are the best objects for setting lim- 
its on violations of the Strong Equivalence Principle 
(SEP) and the Parametrized Post-Newtonian parameter 
«3, which describes Lorentz invariance and momentum 
conservation. These tests use the fact that the gravi- 
tational self-energy of the NS will be much higher than 
that of the WD, and therefore, if the equivalence prin- 
ciples are violated, the two objects will accelerate dif- 
ferently in an external gravitational field or under the 
self-acceleration induced by the velocity relative to a pre- 
ferred reference frame. The net effect on orbits that are 
nearly circular will be to force the ecc entricity into align- 
ment with this acceleration vector l|Damour fc S chafcr 
1991). The prototype of these tests is the se arch for or- 
bital polarization in the Earth-Moon system l)Nordtvedtl 
1968), which currently sets a li mit on the weak-fie l d vi- 
olatio n parameter \n\ of 0.001 l|Dickev et alJll99i IWiil 
2001). The pulsar versions of these experiments test the 
strong-field limit of SEP violation (parameter A) and 
Lorentz invariance/momentum conservation (parameter 
ch), and are thoroughly described in the literature 
lIDamour fc Schaferl Il99lt IWexl 119971 IBell fc Damourl 
H996HWexll200rilStairsll2003USplaver et al.ll2005D . Both 
parameters are identically zero in GR, and 6:3 is predicted 
to be zero by most theories of gravity. 

We now examine the impact of the recently discovered 
binary systems on these tests. The traditional figures of 
merit for choosing systems to test A and &3 are P^/e 
and P b 2 /(Pe), respectively. The other requirements for 
A are that each system must be old enough (ie have 
characteristic age large enough) and must have Co large 
enough that the longitude of periastron can be assumed 
to be randomly oriented; and that each system must have 
Co larger than the rate of Galactic rotation, so that the 
projection of the Galactic accel eration vector onto the or- 
bit can be considered constant ijDamour fc Schaferlll991l: 
!Wexl ll997). Similar requirements hold for d 3 . With its 
extraordinarily low eccentricity, PSR J1853+1303 is a 
prime candidate to help strengthen these tests. The last 
few years have seen the discovery of several other sys- 
tems with compar able or longer-perio d orbits, notably 
PSR,s .72016+1948 (iNavarro et alJ2003l and J0407+1607 
ijLorimer fc Freirdl200a lLorimer et alJl2005|) . We there- 
fore find it worthwhile to update the multi-pulsar anal- 
ysis of Wex (1997, 2000), finding much lo wer limits o n 
each parameter. In keeping with the spirit of lWexl (2000), 
we use all 21 pulsars listed in Table |2 as these are all 
thought to have evolved with similar extended accretion 
periods and therefore represent the overall population of 
such objects. Some of these systems have quite small 
values of P^/e but need to be included nonetheless, as 
possible examples of violation. Our calculation will find 
a median- likelihood value of |A| for each pulsar that cor- 
responds to an induced eccentricity roughly comparable 
to its observed eccentricity, and the combined limit fairly 
represents the limits derivable from the known popula- 
tion. 

We use the following Bayesian analysis. For the SEP 
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A test, we are interested in finding the probability den- 
sity function (pdf) p(\A\ \D,I), where D represents the 
relevant data on the 21 pulsars (namely, their eccentrici- 
ties and longitudes of periastron and associated measure- 
ment errors) and I represents prior information. The un- 
known parameters for each system include the two stellar 
masses, the distance d to the system, and the position an- 
gle f2 of the Line of Nodes on the sky. Given any set of 
these parameters, a "forced" eccentricity vector ep may 
be derived for any given value of A, up to a sign ambigu- 
ity which amounts to flipping the direction of t he vector. 
This can be written <)Damour fc Schaferlll99l : 



|e F | = A 



1 



2FG (mi + m 2 ){2TT / P h f 



(!) 



where c is the speed of light, and, in general relativity, 
F = 1 and G is Newton's constant. Here g± is the 
projection of the Galactic accelerati on vector onto the 
plane of the orbit, and is given by l)Damour fe Schaferl 
Il99l : 



lsx| = | g | [1 — (cost, cos A+sini sin Asin(</>— fi)) ] 



211/2 



(2) 



where <f> is the position angle of the projection of the grav- 
itational acceleration vector g onto the plane of the sky, 
and A is the angle between the line from pulsar to Earth 
and g. Deriving ep requires knowledge of the Galactic 
acceleration at the pul sar position; we assume th e verti- 
cal potential given bv iKuiiken fc Gilmorel l|1989l ) and a 
flat rotation curve with velocity of 222 krns" 1 . A pre- 
diction for the observed eccentricity e b s , P red is then the 
vector sum of ep and a "natural" eccentricity e^. Thus 
the magnitude of eN and the angle 9 between eN and 
— eF are additional parameters. 

For any one j of the 21 pulsars, we may then use Bayes' 
theorem to write: 

p(\A\,i,m 2 ,n i d,e m e\D j ,I) (3) 
<xp(Dj\ \A\,i, m 2 , CI, d, e N , 0, 1)p(\A\,i, m 2 , 0, d, e N , 6\I). 

We compute the integral over eN and 9 separately, ef- 
fectively calculating the marginal p(|A|, i, 7712, ^, d\Dj,I) 
for a particular set of parameters i, m2, and d. For 
each such set of parameters, the angle between ep and 
the true e Q b s , and hence the possible values of eN for 
which the likelihood is significantly non-zero, will almost 
always be very tightly constrained, due to the small mea- 
surement errors on e bs and to. We therefore determine 
the four points representing 3-cr ranges in both e bs and 
lo, and use these to find minimum and maximum values 
of eN and 9. The likelihood for the set of parameters i, 
772,2, and d is then set to 1 for values of eN and 9 fall 
between their minimum and maximum values and oth- 
erwise. Thus the integral will be roughly proportional to 
(6 , ma X -^min)(loge N ,max-fogeN,min), assuming a uniform 
prior on 9 and a log prior on eN- We set the integral to 
zero if eN.min > 0.05 and use 1 x 10~ 6 as a lower bound on 
possible values of eN,min; this conservatively allows each 
pu lsar much more than the eccentricity ranges permitted 
bv lPhinnevI l|1992fl . This approximation to the likelihood 
is necessary as both numerical integration over or (equiv- 
alently) Monte Carlo sampling of the full allowed ranges 
of eN and 9 would be computationally prohibitive. 

For the remaining nuisance parameters i, mi, fi, and 
d we perform a Monte Carlo simulation, drawing m-i 
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Fig. 4. — Posterior probability density functions (pelfs) for the 
test of the Strong Equivalence Principle (SEP). The pdf for each 
pulsar is shown on a linear vertical scale. The horizontal axis is 
displayed logarithmically for clarity, although the range < |A| < 
0.1 was sampled uniformly. The pdfs for some of the pulsars are 
noisy for those cases where the e Q ^ s and gx vectors can be close 
to alignment for certain values of O; these cases are difficult to 
model even with large number of trial systems, but the noise does 
not drastically affect the full pdf. The full pdf p(A\D,I) is the 
normalized product of the individual-pulsar pdfs and is shown on 
a log-log scale in the uppermost right-hand panel. 
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Fig. 5. — Posterior pdfs for the 03 test, similar to the plots in 
Figure [3] 



uniformly from twice the lTauris fc Savoniiei (1999) Pb 
m2 range, cosz from a uniform distribution between 
and 1 (combining these two parameters and the mass 
function yields a value of the pulsar mass mi; we use 
only those systems for which mi is between 1.0 Mq and 
2.5 Mq, or other limits set by timing of the individual 
pulsar), f2 from a uniform distribution between 0° and 
360° and d from a Ga ussian distribution about the best 
iCordes fc Laziol (J2002) distance, assuming a distance un- 
certainty of 25%, or from a Gaussian distribution in par- 
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allax, where measured. For PSR J1713+0747, we restrict 
the parameters i, f2, mi and mi to the region constrained 
by th e recent measureme nt of the orientation of the 
orbit (ISnlaver et alJl20f)l . while for PSR J0437-4715 
we assu me Gaussian distributions about the parameters 
given in Ivan Straten et alJ l)2001l) . We repeat this pro- 
cedure for values of A from 0.00002 to 0.1, in steps of 
0.00002, then normalize; this results in the posterior pdf 
p(|A| \Dj,T) for each pulsar j. For PSR J1713+0747 
alone, the resulting 95% confidence limit o n |A| is about 
0.0158 , similar to the value derived in ISnlaver et al.l 
( 2005). The pulsar data sets are independent, and thus 
we multiply the pdfs to derive pQ A| \D, I). From this, we 
derive a 95% confidence upper limit on|A|of5.6xlO~ 3 . 
Figure0]shows the pdf curves for this test. We note that 
while a logarithmic prior on | A| would result in an upper 
limit a few orders of magnitude smaller, we have chosen 
a uniform prior on |A| in order to be as conservative as 
possible and more consistent with previous work. 

For the 6:3 test, we pro ceed in a similar fashi on. Here 
the forced eccentricity is (|Bell fc Damourll 1996(1 : 



a 3 



^ P |w|P b 2 



24tt P G(mi + m 2 ) 



sin/3 



(4) 



where c p is the gravitational self-energy fraction of 
"compactness" of the pulsar, approximated by 0.21to i 
l)Damour & Esposito-Faresei l992: Be ll fc Damourll 19961. 
and (3 is the (unknown) angle between the pulsar's ab- 
solute velocity w (relative to the reference frame of the 
Cosmic Microwave Background) and its spin vector. For 
this test, we also need the 3-dimensional velocity of the 
system. Where proper motion measurements are avail- 
able, we draw from Gaussian distributions for the proper 
motion to get the transverse velocities; in other cases, 
and always for the unknown radial velocities, we draw 
from Gaussian distributions in each dimension centered 
on the Galactic rotational velocity ve ctor at the pulsa r 
location and with widths of 80 kms -1 l)Lvne et all l998'). 
We sample uniform steps of a% ranging from 1 x 10 2 
to 5 x 10~ 19 . We find a 95% confidence upper limit on 
|d3| of 4.0 x 10 -20 . Figure shows the pdf curves for 
this test. 

The 95% confidence limits wc derive of 5.6 x 10 -3 for 
|A| and 4.0 x 10~ 20 for 1 0:3 1 are considerably better than 
previous lim its of 9 x 10~ 3 and 1.5 x 10~ 19 , respectively 
ll Wexl 12000]) . while still taking into account the contri- 
bution from all pulsars with similar evolutionary histo- 



ries. The SEP test appea rs weaker than the best solar- 
systcm tests of |r?| < 0.001 ijDickev et alll994tlwT!i27)0"l 
but pulsars test the strong-field regime inaccessible to 
the solar-system measurements and are therefore quali- 
tatively different. The | ^3 1 test is nearly thirteen orders 
of magnitude better than tests de rived from the perihe- 
lion shifts of Earth and Mercury l)Willl 11993^ . and again 
tests the strong-field regime of gravity. 

4. CONCLUSIONS 

Of the three new WBMSP systems presented here, at 
least two can be timed at the microsecond level with 
current instrumentation at large telescopes. These ob- 
jects thus show promise for measurement of geometrical 
and/or relativistic timing phenomena in future, giving us 
an idea of the system inclination angles and perhaps the 
masses of the objects. This information will be extremely 
useful in evaluating the validity of the P\r-mi relations 
for estimations of companion masses. All three systems 
should provide proper motion measurements within a few 
years; these will add to our understanding of millisecond 
pulsar velocities throughout their population. Finally, 
in combination with the other low-mass circular-orbit 
systems discovered in recent years, the new pulsars set 
firmer limits on violations of relativistic equivalence prin- 
ciples in the strong-field regime of 5.6 x 10~ 3 for | A| and 
4.0 x 10 -20 for |&s|. A better understanding of the low- 
mass population as a whole will be necessary for further 
improvement of these tests. 
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